EXPONENTIAL RETURN TO EQUILIBRIUM FOR 
HYPOELLIPTIC QUADRATIC SYSTEMS 



M. OTTOBRE, G.A. PAVLIOTIS & K. PRAVDA-STAROV 



Abstract. We study the problem of convergence to equilibrium for evolution 
equations associated to general quadratic operators. Quadratic operators are non- 
selfadjoint differential operators with complex-valued quadratic symbols. Under 
appropriate assumptions, a complete description of the spectrum of such opera- 
tors is given and the exponential return to equilibrium with sharp estimates on 
the rate of convergence is proven. Some applications to the study of chains of 
oscillators and the generalized Langevin equation are given. 



1. Introduction 

The spectral properties of hypoelliptic operators and the problem of return to 
equilibrium for the associated parabolic equations are topics of current interest. Such 
operators arise naturally in non-equilibrium statistical mechanics as the generators of 
Markov processes which model the dynamics of classical open systems [3J [5J [HI [551 US] • 
A standard example of hypoelliptic equation is the Fokker-Planck equation [30] , 

(1) f d t f + v-V x f~V x V{x)-V v f- 1 V v -(v + P~ 1 V v )f = {), 
\ f(x,v,t)\ t=0 = f {x,v). 

This equation is the forward Kolmogorov equation for the Langevin stochastic differ- 
ential equation 

(2) ( x = v, 

where W(t) denotes a standard Brownian motion in R d , 7 the friction coefficient and B 
the inverse temperature. Several properties of the Fokker-Planck equation have been 
studied including spectral properties and convergence to equilibrium [3J[H1[II], homog- 
enization H] , short time asymptotics [TUJ and semiclassical limits [TT], [T^l E3J E5J . 
Related problems have been studied for the Fokker-Planck equation corresponding to 
the generalized Langevin equation under a Markovian approximation [25] . 

Recently, the theory of hypocoercivity [33] has been developed for studying the 
convergence to equilibrium of kinetic equations. This theory applies to evolution 
equations of the form 

(3) d t f + (A*A + X )f = 0, 

where A, X$ are some operators defined on an appropriate Hilbert space H, which 
is usually the L 2 space weighted by the invariant measure of the underlying Markov 
process, with X being skew-adjoint in this space, A"q = — X Q . Under appropri- 
ate assumptions on the commutators between the operators A, Xq and higher order 
commutators, it is possible to prove the exponential return to equilibrium for these 
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evolution equations. Regarding the Fokker-Planck operator, quantitative estimates 
on the rate of convergence can be obtained via the analysis of the spectral gap. How- 
ever, it is not possible in general to determine exactly the spectrum and one of the 
few instances when this is possible is the case of quadratic potentials [T5] . 

In the present work, we consider evolution equations associated with general qua- 
dratic operators 

(4) | ^(t,x) + q w (x,D x )u(t,x) = 

I u(t,-)\ t=0 = u Q eL 2 {K n ), 

and address the problem of the exponential return to equilibrium for these systems. 
Quadratic operators are pseudodifferential operators defined in the Weyl quantization 

(5) q w {x,D x )u{x) = 1 ^ n [ e«*-Mq(^-,t)u{y)dydt, 

by symbols q(x,£), with (x, £) G R" x K" and n G N*, which are complex- valued 
quadratic forms 

ij : R" x R" ^ C. 

These operators are differential operators with simple and fully explicit expression. 
Indeed, the Weyl quantization of the quadratic symbol x a ^, with (a, /J) G N 2 ™, 
| a + ;9| = 2, is the differential operator 

^ — x - — , D x = i 1 d x . 

Notice that quadratic operators are non-self adjoint operators and that those with 
symbols having non-negative real parts are accretive. Since the classical work by 
J. Sjostrand [3T], a complete description for the spectrum of elliptic quadratic opera- 
tors is known and has played an important role in the analysis of partial differential 
operators with double characteristics. Elliptic quadratic operators are quadratic op- 
erators whose symbols satisfy the ellipticity condition 

(6) (x,OeM 2 ", q(x,O=0^(x,O = 0. 

In a recent work |16j . the spectral properties of non-elliptic quadratic operators, that 
is operators whose symbols may fail to satisfy the ellipticity condition on the whole 
phase space M 2n , were investigated. For any quadratic operator whose symbol has a 
real part with a sign, say here a symbol with non-negative real part 

(7) Re q > 0, 

it was pointed out the existence of a particular linear subvector space S in the phase 
space E 2 ", S C M 2 ", intrinsically associated to the symbol q and called singular space, 
which plays a basic role in the understanding of the properties of this non-elliptic 
quadratic operator. In particular, this work |16j (Theorem 1.2.2) gives a complete 
description for the spectrum of any non-elliptic quadratic operator q w (x,D x ) whose 
symbol q has a non-negative real part, Re q > 0; and satisfies an assumption of partial 
ellipticity along its singular space S, 

(8) (i l OeS )9 (i ) fl=fl=>(i,0=0. 

Under those assumptions, the spectrum of the quadratic operator q w (x, D x ) is showed 
to be composed of a countable number of eigenvalues with finite multiplicity and the 
structure of the spectrum is similar to the one known for elliptic quadratic opera- 
tors [31] . This condition of partial ellipticity is weaker than the condition of ellipticity 
in general, S C M. 2n , and allows to deal with more degenerate situations. An impor- 
tant class of quadratic operators satisfying condition (Tsl) are those with zero singular 
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space S = {0}. In this case, the condition of partial ellipticity trivially holds and 
these differential operators have been showed to be hypoelliptic [27] (Theorem 1.2.1). 
More precisely, they enjoy specific subelliptic properties with a loss of derivatives 
with respect to the elliptic case depending directly on the structure of their singular 
spaces. We refer the reader to [27JI2E] for a complete presentation of these subelliptic 
properties and examples of such subelliptic quadratic operators. 

In the present work, we shall consider this class of hypoelliptic operators with zero 
singular space and study further the structure of the bottom of their spectra. More 
specifically, we shall see that the first eigenvalue in the bottom of their spectra has 
always algebraic multiplicity one with an eigenspace spanned by a ground state of 
exponential type. We shall also give an explicit formula for the spectral gap which 
is computable via a simple algebraic calculation; and finally answer the question of 
long time behavior of the associated evolution equations by proving the property of 
exponential return to equilibrium for these quadratic systems. 

Let us finally mention that the singular space theory [THl [TTJ [TSl [STJ apply in 
various settings and that certain applications to the study of chains of oscillators and 
the generalized Langevin equation are discussed at the end of this article. 

2. Statements of the main results 

We begin by recalling miscellaneous facts and notations about quadratic operators. 
In all the following, 

g:l"xlj -> C 

(x,t) ' ^ q(x,0, 
stands for a complex-valued quadratic form with a non-negative real part 

(9) Re q(x,£) > 0, (af,£) € E 2n ,n G N*. 

We know from [5D] (p. 425-426) that the maximal closed realization of the operator 
q w (x,D x ), i.e., the operator on L 2 (R n ) with domain 

D(q) = {u G L 2 (R n ) : q w {x,D x )u G L 2 (R n )}, 

coincides with the graph closure of its restriction to S(M. n ), 

q w (x,D x ) : S(R n ) -> S(R n ); 

and that this operator is maximally accretive 

Re(q w (x,D x )u,u) L2 > 0, ueD(q). 

Associated to the quadratic symbol q is the numerical range £(<?) defined as the 
closure in the complex plane of all its values, 

(10) E(g) = g(ffi£ x KJ?). 

The Hamilton map F G Mi n (C) associated to the quadratic form q is the unique map 
defined by the identity 

(11) q((x,0;(y,v))=<r((x,0,F(y,ri)), Me^'M^eK 2 "; 

where <?(•; •) stands for the polarized form associated to the quadratic form q and a 
is the canonical symplectic form on R 2n , 

(12) a((x, 0, (y, rj)) = £.y - x.r,, (x, G R 2n , (y, V ) G R 2n . 
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It directly follows from the definition of the Hamilton map F that its real part and 
its imaginary part 

Re F = ~(F + F) and Im F = —(F-F), 

are the Hamilton maps associated to the quadratic forms Re q and Im q, respectively. 
One can notice from ( 11 ) that a Hamilton map is always skew-symmetric with respect 
to a . This is just a consequence of the properties of skew-symmetry of the symplectic 
form and symmetry of the polarized form 

(13) VX, Y G E 2 ™, a{X, FY) = q(X; Y) = q{Y; X) = a(Y, FX) = -a{FX, Y). 

Associated to the quadratic symbol q is the singular space S defined in [16] as the 
subvector space in the phase space equal to the following intersection of kernels 

(14) S= ( p| Ker[Re F(Im F) J ]) fli 2 ", 

3=0 

where the notations Re F and Im F stand respectively for the real and imaginary 
parts of the Hamilton map associated to q. Notice that the Cayley-Hamilton theorem 
applied to Im F shows 

(Im F) k X G Vect(X, (Im F) 2n ~ 1 X) , X G R 2n , keN; 

where VectfX, (Im F) 2n ~ 1 X) is the vector space spanned by the vectors X,..., 
(Im F) 2n ~ 1 X. The singular space is therefore equal to the finite intersection of the 
kernels 

2n-l 

2n 



(15) S = ( P| Ker[Re F(Im F) j ] 



n. 



As mentioned above, when the quadratic symbol q satisfies a condition of partial 
ellipticity along its singular space S, 

(16) (i,Oe S, q(x,t) =0=>(as,O=0, 

Theorem 1.2.2 in [16] gives a complete description for the spectrum of the quadratic 
operator q w (x, D x ) which is only composed of eigenvalues with finite algebraic multi- 
plicity 

(17) a(q w (x,D x )) = { ^ (r A + 2k x ) (-i\) : k x G n}, 

-iAec+u(s( 9 | s )\{o}) 

where F is the Hamilton map associated to the quadratic form q, r\ is the dimension 
of the space of generalized eigenvectors of F in C 2n belonging to the eigenvalue A G C, 

T,(q\ s ) = qJS) and C + = {zeC:Rez>0}. 

In the present paper, we study evolution problems associated to accretive quadratic 
operators q w (x, D x ) with zero singular space S = {0}. Notice that in this case q(S) = 
{0} and the condition of partial ellipticity along the singular space trivially holds. 
The spectrum is then reduces to 

(18) a(q w (x,D x )) = { Yl (rx + 2k x )(-zX):k x eNy 

Xea(F) 
-iAeC + 
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Define 

(19) Mo= - lXr * 

Xea(F) 

-i\ec + 

and 

(20) t = min Re(2(-tA)) = 2 min Im A, 

Ae<r(F) V ' \£a(F) 

-iAeC + Im A>0 

Our first result shows that the first eigenvalue in the bottom of the spectrum /iq 
has always algebraic multiplicity one with an eigenspace spanned by a ground state 
of exponential type and that the spectral gap is exactly given by the positive rate 
r > 0. 

Theorem 2.1. Let q w (x,D x ) be a quadratic operator whose symbol 

q : R£ x R£ C 

(x,£) H> q{x,Cl, 

is a complex-valued quadratic form with a non-negative real part Re q > and a zero 
singular space S = {0}. Then, the first eigenvalue in the bottom of the spectrum 



defined in ( 19 I has algebraic multiplicity one and the eigenspace 

Ker(q w (x,D x )- IM) )=Cu , 

is spanned by a ground state of exponential type 

u (x) = e~ a ^ G 5(M n ), 

where a is a complex-valued quadratic form on M. n whose real part is positive definite. 
Moreover, the spectral gap of the operator q w (x,D x ) is exactly given by the positive 
rate t > defined in p0| , 

a(q w {x, D x j) \ {fi Q } C {z G C : Re z > Re ^ + r }. 



Let us underline that operators satisfying the assumptions of Theorem |2.1| are not 
elliptic in general and that their symbol may fail to satisfy the ellipticity condition 

(x,OeR 2n , q(x,0=0^(x,0=0. 

However, the assumption about the singular space S = {0} ensures that such an 
operator enjoys noticeable subelliptic (hypoelliptic) properties 

(21) || {((x,0) 2{1 - 5 T4l^ < C(\\q w (x,D x )u\\ L2 + \\u\\ L2 ), 

with ((x, £)) = (1 + \x\ 2 + ICI 2 ) 1 ^ 2 and a loss of derivatives with respect to the elliptic 
case < 5 < 1. More specifically, consider 

q :R X x M.£ -> C 

(x,Z) ^ q{x,0, 

a complex-valued quadratic form with a non-negative real part 

Re q(x,£) > 0, (x,£) G M 2 ", n G N*, 

and zero singular space S — {0}. According to the definition of the singular space 



( 15 ), there exists a smallest integer < ko < 2n — 1 such that 

k 

(22) ( p| Ker[Re F(Im F) J ] J n R 2n = {0}, 

3=0 
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FIGURE 1 . Notice that the first eigenvalue in the bottom of the spec- 
trum is not necessarily real. 

Im zl 



T 



»/'« 



Rc z 



where F stands for the Hamilton map of q. Theorem 1.2.1 in [57] shows that the 
quadratic operator fulfills a global subelliptic estimate with loss of 8 = 2ko/(2ko + 1) 
derivatives with respect to the elliptic case, that is, there exists C > such that for 
all u G D(q), 

(23) ||(<(s,0> 2/(a *° +1) r«IL a ^ C(\\q w (x,D x )u\\L* + \\u\\v), 

with 

D(q) = {u £ L 2 (R n ) : q w {x, D x )u G L 2 (R n )}. 

Those subelliptic estimates are a key argument for getting sharp resolvent estimates 
and proving the property of exponential return to equilibrium. Next result proves 
this property of exponential return to equilibrium and provides an exact formula for 
the optimal rate of convergence: 



Theorem 2.2. Let q w (x,D x ) be a quadratic operator satisfying the assumptions of 
Theorem 2.1 By using the notations introduced in (19) and (20), we consider the 
operator 

Q = q w {x 1 D x ) - hq. 
Then, for all < r < To, there exists a positive constant C > such that 

Vt>0, ||e- tc3 -n |U( i2 ) <Ce~ rt , 

where Ho is the rank-one spectral projection associated with the simple eigenvalue zero 
of the operator Q and || • ||,c(l 2 ) stands for the norm of bounded operators on L 2 (M. n ). 
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Let q w (x,D x ) be a quadratic operator satisfying to the assumptions of Theo- 
rems |2.1| or |2.2[ and /io be the first eigenvalue in the bottom of the spectrum. Theo- 
rem [2J] shows that 

A^o = -«Ar A G C, 

AGcr(_F) 
-iAeC+ 

where F stands for the Hamilton map of the quadratic form q. We shall now consider 
the specific case when this quadratic operator is real in the sense that q w (x, D x )u is a 
real-valued function whenever u is a real-valued function. We begin by checking that 
this assumption ensures that /io is real and that the quadratic form a defining the 
ground state 

u (x) = e~ a ^ G S(R n ), q w (x, D x )u = » u , 
is then positive definite. Indeed, by using that the operator is real and passing to the 
complex conjugate in the equation 

q w (x,D x )e- a W =fi e- a M, 

we obtain that 

q w (x,D x )e-"W =Ji e-< x \ 

It shows that /Z is also an eigenvalue for the operator q w (x,D x ). Since by Theo- 
rem [2T] we have 

a(q w (x,D x )) \ {/i } C {z G C : Re z > Re /i + r }, 

with t > 0, Re /i = Re ~p implies that the eigenvalue /io is real. Furthermore, 
notice that e~ a and e~ a are therefore two eigenvectors associated to the eigenvalue 
/io- We deduce from Theorem 2T that there exists Ai, A2 G C, not both equal to zero 
such that for any x G E™, 

Recalling that a is a quadratic form, we obtain with a; = that Ai = — A2, which 
implies successively that e~ a — e~ a = 0, e~ lb — e lb — 0, e 2lb — 1, sin(26) = 0, 
cos(26) = 1, if b stands for the imaginary part of the quadratic form a. This shows 
that the quadratic form b is constant 

= d Xj ( sin(26)) = 2 cos{2b)d x . b = 2d x . 6. 

and therefore identically equal to zero. The quadratic form a is therefore a positive 
definite quadratic form. 

On the other hand, notice that the adjoint operator (see [50), p. 426) is actually 
given by the quadratic operator 

q w (x,D x y =q w (x,D x ), 

whose symbol is the complex conjugate symbol of q. It is therefore a complex-valued 
quadratic form with non-negative real part and a zero singular space for which Thc- 



2.1 applies. 



Then, /io is also the first eigenvalue in the bottom of the spectrum for the qua- 
dratic operator q w (x, D x )* , and we know from Theorem 2.1 that the eigenspaces 
associated with the eigenvalue /i for both operators q w (x,D x ) and q w (x,D x )* are 
one-dimensional subvector spaces with ground states of exponential type. We shall 
assume further that the two operators have same ground state 

(24) Kei(q w (x, D x ) - Mo ) = Ker(<f>, D x )* - Mo ) = Cu Q C S(R n ), 

with uq(x) — e~ a ( x \ x £ R n ; where a is a positive definite quadratic form on K™. 



Under those assumptions, one can specify further the result given in Theorem 2.2 
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Theorem 2.3. Let q w (x,D x ) be a quadratic operator satisfying the assumptions of 



Theorem 2.1 Assume that this quadratic operator is real and satisfies (24). Using 



the notations introduced in (19 1 and (20), we consider the operator 

Q = q w (x,D x )~fl . 
Then, for all < r < To, there exists a positive constant C > such that 



Vt > 0,VuG L 2 (R n ), \\e- ttJ u-c u u \\ L 2 {Rn) < Ce' Tt \\u\\ L 2 {mnh 

2 

L 2 (R™)> 

2 



where c u is the L 2 (R n ) scalar product of u and wo/||W 2 



°u = ||Mo|| L 2( R „)(w, uo)l 2 (R")- 

We finally recall that the singular space theory JTHl HZ1 UH EH1 HE] was mainly 
inspired by the seminal work on the Fokkcr-Planck equation |15j and aims at extending 
to wider classes of doubly characteristic operators some results proved for Fokker- 
Planck operators in [TTJ [TU EH E3- After preliminary advances in this direction 
in [2B], the main algebraic fact discovered in [TH] was to notice that the symbols of 



operators satisfying the assumptions of Theorem 2.1 have real parts which may not 
be positive definite in general; however, when the real part of these quadratic symbols 
are averaged by the flow of the Hamilton vector field associated to their imaginary 
parts Hi mq , 

(25) (Re q) T {X) = -L J Re q(e tH ^X)dt » 0, 

with X = (x,£) £ M. 2n , quadratic forms (Re q)x become positive definite for any 



T > 0. This averaging condition (25 ) is a key property in the works (TTJ rT2J [T31 ITS] . 
One of the main interest of the singular space theory is to provide an easily computable 
algebraic condition on symbols of general quadratic operators for checking that the 



averaging property (25) holds. The remainder of this article is organized as follows. 
Next section is devoted to the proofs of Theorem |2.1[ |2.2| and |2.3| whereas various 
applications of these results to the study of chains of oscillators and the generalized 
Langevin equation are discussed in Section [4] 

Acknowledgement. The third author would like to express his gratefulness to 
Johannes Sjostrand for enriching remarks and comments on the present work. 

3. Proof of the main results 

This section is devoted to the proofs of Theorems |2.1| |2.2| and |2.3| 

3.1. Proof of Theorem |2.1[ We begin by noticing that the set 

(26) {A £ C : A e a(F), -iX G C+} = {A € C : A € a(F), Im A > 0}, 

is non-empty. Indeed, the work |16) gives an alternative description of the singular 
space in terms of eigenspaces of the Hamilton map F associated to its real eigenvalues. 
More specifically, when q is a complex-valued quadratic form with a non-negative real 
part fulfilling the partial cllipticity condition on its singular space 

(s,0eS> q(x,t)=0=>(x,t)=0, 

which is obviously satisfied when S = {0}, then the set of real eigenvalues of the 
Hamilton map F can be written as 

cr(F)nR={A 1 ,...,A r ,-A 1 ,...,-A r }, 
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with \j 7^ and Xj ^ ±Xk if j ^ k. The singular space is then ([Mj, p. 807) the direct 
sum of the symplectically orthogonal spaces 

(27) S = S Xl (B a± S X2 O a± ...® a± S Xr , 
where S\ } is the symplectic space 

(28) 5 Aj = (Kcr(F - Xj) Ker(F + Xj)) n R 2n . 
Recalling from Proposition 4.4 in [50] that the subvector space 

Ker(F — Xj) ® Ker(F + Xj), 

is the complexification of its intersection with M 2 ™, condition S = {0} therefore implies 
that the Hamilton map F has no real eigenvalue 

a(F) fll = 0. 

On the other hand, we know from [T5] (Lemma 21.5.2) that V\ the space of generalized 
eigenvectors of F belonging to the eigenvalue A e C is the dual with respect to 
the symplectic form of the space of generalized eigenvectors V-\ belonging to the 
eigenvalue —A. It follows that 

A G a{F) O-Ae <t(F). 

This implies that the set 

{A G C : A G a(F), -iX G C+} = {A G C : A G a(F), Im A > 0}, 



is non-empty; and that the quantities /Uo G C and tq > are well-defined in ( 19 ) and 
(20 1. Since from (18) the spectrum of the operator q w (x,D x ) is only composed of the 
following eigenvalues with finite algebraic multiplicity 

(29) a(q w (x,D x )) = { £ (r A + 2fc A ) (-iX) : fe A € n}, 

AGcr(_F) 
-i\eC + 

we notice that the first eigenvalue in the bottom of the spectrum is equal to 

-i\ec + 

and we easily identify the spectral gap to be the positive rate 

t = min Re(2(— iXj) = 2 min Im A > 0. 

AG<t(F) V ' \£a(F) 

-i\£C+ Im A>0 

It remains to check that this eigenvalue /io has algebraic multiplicity one with an 
eigenspace 

Ker(q w (x,D x )-M) =Cu , 
spanned by a ground state of exponential type 

U (x) = e~ a ^ G S(R n ), 

where a is a complex- valued quadratic form on M." whose real part is positive definite. 

In order to do so, we shall first consider the case when the real part of the quadratic 
symbol Re q is positive definite. In this specific case, the symbol q is elliptic on the 
whole phase space M 2 ™, 

(x,0eR 2 ", q(x,O = 0^(x,O = 0; 
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and one can refer to the original work by J. Sjostrand on elliptic quadratic differential 
operators [31] . Following this fundamental work on which relies the analysis led in [16] , 
we consider the positive Lagrangian plane (Proposition 3.3 in [31)). 

A£<r(F) 
-iAeC+ 

where V\ is the space of generalized eigenvectors belonging to the eigenvalue A of the 
Hamilton map F associated to the quadratic symbol q. We recall that a complex 
subspacc A C C 2 " of dimension n is called Lagrangian if the canonical symplectic 
form a vanishes identically on A. Furthermore, a Lagrangian plane A C C 2n is said 
to be positive if 

-ia(X,X) > 0, 

for all 7^ X = (x,tf) 6 A. According to (3.5) in [3T] (p. 95), there exists a symmetric 
matrix B + whose imaginary part Im B + is positive definite such that 

V+ = {{x,B+x) : x e C"}. 

Define the two quadratic forms 

b + (x) = (x,B+x), x e R n ; a(x) = ~ib + (x), x e R". 

Notice that the real part of the quadratic form a is positive definite. Keeping on 
following the analysis led in [31] (p. 98), on can then use a Jordan decomposition for 
the Hamilton map F in order to find a canonical transformation IC of C 2 " such that 

rx r\- l 

$ = 2 H {Yl Xx iMj,* + ' -\) , 

\£<t(F) j=l 3=1 
-i\eC + 

with q = qoK" 1 . Here, Xj,\,£j,\, with 1 < j < r\, A € cr{F) and — iX G C+; stands for 
a relabeling of the standard real symplectic coordinates (x, £) € M 2 ™; and -fj.\ £ {0, 1} 
with 7o,a = 7r A ,A = 0. We recall that r\ stands for the dimension of V\ the space of 
generalized eigenvectors associated to the eigenvalue A G C. This reduction is proved 
in [3T] (p. 98) in (3.15) and (3.16). Consider the two quadratic operators 

\ea{F) 3 = 1 

-i\ec+ 

and 

q£(x,D x )=2 X! ~,>..v0.a/- ) : . v 

AGcr(F) 3 = 1 

Notice that these two operators commute and that 

q w (x,D x )=q™(x,D x ) + q™(x,D x ). 

Proposition 3.10 in [31] shows that every monomial is a generalized eigenvector of the 
operator 

q w (x,D x );P(R n )^P(R n ), 

where P(R n ) stands for the space of complex polynomials regarded as functions on 
R. n . More specifically, the spectrum of the operator q w (x,D x ) acting on the space 
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P(R n ) is exactly given 

(30) { J2 (rA+2Afc)(-iA):Afc€N}. 

\ecr(F) 

-i\ec + 

Proposition 3.10 in [3T] also shows that the space of generalized eigenvectors 

E (q w (x, D x ) , P(R n ), fi) , 
associated to an eigenvalue /x is of finite dimension 

dim E(q w (x,D x ),P(R n ),fi) < +oo. 
Inspecting the proof of Proposition 3.10 in [31] shows that 

2J \(x jt \D x . x + D x ■'' ;..\ j /'.\ = -«A(r A + 2fc A )pA, 

3=1 

for any homogeneous polynomial of degree k\ in the variables Xi t \,..., x r (\\\. This 
implies that the monomials 

P(x) = Yi P*( x )' 

AGo-(F) 

-z\ec + 

constitute a basis of eigenvectors for the operator q™(x,D x ) acting on P(R"), 
q?(x,D x )p= {rx + 2k x )(-i\)p. 

Aeo-(F) 
-iAeC+ 

According to this description, we notice that the first eigenvalue in the bottom of the 
spectrum for the operator q^(x, D x ) acting on the space P(K") is the value /io defined 
in (19). Furthermore, this eigenvalue /io has algebraic multiplicity 1. Indeed, all the 
non-negative integers k\ must necessarily be zero and the only possible eigenvector 
associated to /iq is up to a non-zero constant the constant polynomial 1. Keeping 
on following the proof of Proposition 3.10 in [31j . one may notice that the operator 
q2( x iD x ) is nilpotcnt on any space of generalized eigenvectors 

E(q?(x,D x ),P(pL n ) t n), 

with (USC; for the operator q™(x, D x ) seen as acting on the space P(M n ). This proves 
in particular that /iq is an eigenvalue with algebraic multiplicity 1 for the operator 

q w (x,D x ):P(« n )^ P(R n ), 

with eigenvector equal to the constant polynomial 1. According to [31] (p. 100), there 
exists a bijection 

K : P(R n )e ib+{x)/ ' 2 -> P(M"), 

such that 

(31) q w (x,D x )oK = Koq w (x,D x ). 

More specifically, according to (3.22), (3.23), (3.24) and (3.25) in [31 , this mapping 
may be written as 

(32) K = K 3 K 2 K X : P(R n )e lb+ -> P(R n ), 
with some operators 

(33) Ki : P(R n )e lb+{x)/2 -> P(M n ), K lU = e- lb+{x)/2 u; 

(34) K 2 : P(R n ) -> P(R n ), K 2 u = F' 1 (e ib ~ Fu) ; 
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and 

(35) K a ■ P(M. n ) -> P(R n ), (K 3 u)(x) = u{Cx); 

where J- is the Fourier transformation, b~ a quadratic form and C an invertible 
complex matrix. Proposition 3.11 in |31j then establishes that the spectrum of the 
operator 

q w (x,D x ) : P(R B )e a+ W/ 2 -> P(R n )e i6+ ^/ 2 , 
is also exactly given by (30); and that for every eigenvalue the corresponding space 
of generalized eigenfunctions is finite dimensional. More specifically, we deduce from 
(31), (32), (33), (34) and (35) that /Jo is an eigenvalue with algebraic multiplicity 1 
for the operator 

q w (x,D x ) : P(R n )e lb+{x)/2 -> P{R n )e lb+{x)/2 , 
with associated eigenvector K~~ 1 {1). An explicit computation shows that 

with a = —ib + /2. As a final step in the proof of Theorem 3.5, it is proved in [3T] 
(p. 102) that for any \i £ C, the space of generalized eigenvectors associated to the 
eigenvalue /i for the operator 

q w (x,D x ) : P(R n )e lb+{x)/2 -> P(R n )e lb+{x]/2 , 

exactly coincides with the space of generalized eigenvectors associated to the eigen- 
value fj, for the operator 

q w (x,D x ) :L 2 (R n )^L 2 (R n ), 

that is 

E(q w (x, D x ),P(R n )e lb+ / 2 , /*) = £(<T(x, D x ), L 2 (R n ), M ) . 
This implies that when the real part of the quadratic symbol Re q is positive definite, 
Ho the first eigenvalue in the bottom of the spectrum of the operator q w (x, D x ) acting 
on L 2 (R n ) has algebraic multiplicity 1 with a ground state of exponential type 

q w (x,D x )u = nou , 

with Uo(x) = e~ a ( x ^ € <S(K n ), where a is a complex-valued quadratic form on R n , 
whose real part is positive definite. This proves Theorem 2.1 in the special case when 
the real part of the quadratic symbol Re q is positive definite. 
Let us now consider the general case when the quadratic symbol 

q : R^ x R£ C 

0,0 ^ q(x,0: 

is a complex-valued quadratic form with a non-negative real part Re q > and a zero 
singular space S = {0}. We shall now inspect the analysis led in [16] in order to prove 
Theorem |2.1| in the general case. Notice that in this general case the real part of the 
quadratic symbol Re q is not necessarily positive definite. However, as noticed in [IB] 
(p. 809-810) (see also Proposition 2.0.1 in [T^]) and mentioned above, the fact that the 
singular space S is reduced to zero implies that the average of the real part of the 
quadratic form q by the flow generated by the Hamilton vector field of its imaginary 
part Hi mq , 

(Re q) T {X) = -L J ^ Re q (e tH ^X)dt, X = (x,0 € R 2n ; 

is positive definite for any T > 0. We shall now inspect more specifically the proof 
of Proposition 3.1.1 in jl6j. Notice that we consider in [16] complex- valued quadratic 
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symbols with non-positive real parts. One can therefore apply all this analysis when 
considering the symbol — q instead of q and this accounts for sign differences when 
referring to results proved in [IB]. Following [IB] (p. 826), we shall discuss certain 
IR-deformations of the real phase space M. 2n , where the averaging procedure along 
the flow defined by the Hamilton vector field of Im q plays an important role. To that 
end, let G = Gt be a real-valued quadratic form on R 2n satisfying 

(36) H lmg G = -Req+(Req) T . 



As in [TT] , we solve ( 36 ) by setting 



(37) G(X) = J k T (t)Re q{e tHl ™* X)dt, 

where k T (t) = k(t/2T) and k e C(R \ {0}) is the odd function given by 
k{t) = for |t| > - and k'(t) = -1 for < \t\ < 

Let us notice that k and kx have a jump of size 1 at the origin. Associated with G is 
a linear IR-manifold, defined for < e < eo, with e > small enough 

(38) A eG = e l£liG (R 2n ) c C 2n , 

where e leH ° stands for the flow generated by the linear Hamilton vector field ieFIq 
taken at time 1. Keeping on following I16| (p. 827) (see (3.1.23)), we consider the 
quadratic symbol 

(39) q s {X) = q(e ieH °X), 
whose real part is easily seen to be equal to 

Re q £ (X) = Re q(X) + eH lmq G(X) + 0(e 2 \X\ 2 ). 



We deduce from (36) and the facts that the quadratic form (Re q)r is positive definite 
and Re q > 0, 

(40) Re q e (X) = (1 - e)Re q(X) + e(Re q) T {X) + Oie 2 \X\ 2 ) > J \X\ 2 , 

for C > 1 and 0<e<l. We proved in jTEj (Lemma 3.1.1) that for any < e < eq, 
with Eq > sufficiently small, the spectrum of the quadratic operator q™(x,D x ) is 
actually independent of the parameter < e < Sq, 



V0 



<e<e , a(q™{x,D x )) = { ]T (r x + 2k x )(-i\):k x €®}, 



\<£t{F) 



where F is the Hamilton map associated to the quadratic form q, r\ is the dimension of 
the space of generalized eigenvectors of F in C 2 ™ belonging to the eigenvalue A G C. 
The independence of the spectrum of the operator q™(x,D x ) with respect to the 
parameter < e < £o, is linked to the fact that the Hamilton map F £ of the quadratic 
form q e is isospectral for any < e < 1 with the Hamilton map F of the quadratic 
form q, since the symbols q and q e are related by a canonical transformation. As 
a final step in [TB], we proved in (p. 830-831) that if the positive constant eo > is 
chosen sufficiently small then the spectrum of the operator q w (x,D x ) is equal to the 
spectrum of the operator q™(x,D x ) for any < e < eo- Furthermore, eigenvectors 
and generalized eigenvectors agree, see [TB] (p. 830-831) and [T5] (Proposition 2.1). 
This implies that one can deduce the result of Theorem |2.1| in the general case when 
the quadratic symbol has a non-negative real part Re q > and a zero singular space 
S = {0} from the one proved previously in the particular case when the real part of 



the symbol is positive definite. This ends the proof of Theorem 2.1 
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3.2. Proof of Theorem 12.21 Wc consider 

q : R£ x RJ? C 

a complex-valued quadratic form with a non-negative real part 

Re q(x,£) > 0, G K 2n , nef, 

and zero singular space S 1 = {0}. According to the definition of the singular space 



( 15 ), there exists a smallest integer < fco < 2n — 1 such that 

fe 

(41) I p| Ker[Re F(Im Fy]j n M 2 " = {0}, 

j'=o 

where F stands for the Hamilton map of q; and Theorem 1.2.1 in [57J shows that the 
operator q w (x,D x ) fulfills a global subelliptic estimate with loss of S = 2fco/(2fco + 1) 
derivatives with respect to the elliptic case, that is, there exists C > such that for 
all u G D(q), 

\\(((x,0) 2/{2ko+1 T u \\^ ^ C(\\q w (x,D x )u\\ L z + \\u\\ L2 ), 

with 

D(q) = {ue L 2 (R n ) : q w (x, D x )u G L 2 {R n )}. 
By using the proof given for establishing Theorem 1.2.1 in [57], one can actually 
directly obtain the more general estimate that there exists a constant C > such 
that for all u G D{q) and v G K, 

(42) || 0) 2/i2ko+1 T4 L i ^ C(\\q"(x, D x )u - ivu\\ L , + \\u\\ L ,). 

To prove this fact, it is actually sufficient to notice that one can replace the operator 
gWick k v the operator g Wlck — iv in the estimate (2.10) in [27]. Indeed, we recall]^] 
that real Hamiltonians get quantized in the Wick quantization by formally selfadjoint 
operators on L 2 . Since the weight function g given by Proposition 2.0.1 in |27j is 
real-valued and s a positive constant, we have the equality 

Re((l - eg) Wick (q Wick - iv)) = Re((l - e 5 ) Wick g Wick ) , 

which justifies that one can actually replace the operator g Wlck by the operator g Wlck — 
iv in the estimate (2.10) in [27] . Then, by noticing that 

Re{q wick -tv) =Re((z Wick ), 

one can also replace the operator g Wlck by the operator q Wlck — iv in the estimate 
preceding the estimate (2.19) in [57]. Finally, after having done these two slight 
modifications, one can then use exactly the same proof as the one given for proving 



Theorem 1.2.1 in [27] in order to establish the estimate (42). 

We now need to recall few facts about the relationship between pseudo-differential 
calculus and functional calculus when using the Sobolev scale A r , r G K; defined by 
the operator 

A 2 = (i + |x| 2 + iei 2 r-i + ix| 2 + i^i 2 . 

In order to do so, it will be convenient to directly refer to the presentation given by 
F. Herau and F. Nier in [TJ] (Appendix A, p. 205). Our study corresponds to the 
specific case when the potential V in [TJ] is taken quadratic 

V{x) = \\x\ 2 . 



1. see the appendix about Wick calculus in [27] 
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Define 



and 



{u e S>(R n ) ■■ (((x,0) r ) w u e L 2 (R n )}, 



r e 



ll«lk = ||«( i r,0) r r«ll^. 
By referring to [Hj (Proposition A. 4), we notice that, for any r £ M, the operator 
A r is a pseudo-differential operator whose Weyl symbol belongs to the symbol class 
S(((x^)} r ,dx 2 + d£ 2 ) composed by C°°(R X ^,C) functions satisfying 

Va e N 2n ,3C a > 0,V(a:,O £ R 2n , \d^a{x,^)\ < C a ({x,OY- 

Furthermore, for any r £ R, the domain of the operator A r is equal to the space H r 
with equivalence of the two norms and ||A r u||i2. It follows from (42 1 that there 

exists a constant C > such that for all u £ D(q) and ygl, 

(43) \\A 2 ^ 2k ^u\\ 2 L2 < C(\\q w (x, D x )u - iuu\\] 

Consider now the metric 

dx 2 + d£ 2 



\\h)- 



and weight M(x,f) = ({x,0) 2 



We recall for instance from [22] (Lemma 2.2.18 



and section 2.2) that this metric is admissible with gain M. Any quadratic form is 
obviously a first order symbol in the symbolic calculus associated to this metric 

qeS(M,g). 

One can then directly deduce from the Fcfferman-Phong inequality (See for instance 
Theorem 2.5.5 in [55]) that there exists a positive C > such that for all u £ D(A 2 ), 

(44) \\q w (x,D x )u\\ 2 L2 <C\\A 2 u\\ 2 L2} 

with 

D(A 2 ) = {u£ L 2 (R n ) : (1 + \x\ 2 + \D x \ 2 )u £ L 2 (R n )}. 
Notice from Theorem 2.1 together with ( p~9[ ) and (26) that 

a(q w {x, D x )) n iR = a disc (q w (x, D x )) n « = 0. 

One can apply the abstract functional analysis led in [S] Sec. 6.1] in order to obtain a 
specific control of the resolvent of the operator q w (x, D x ) in particular regions of the 
resolvent set. More specifically, following the analysis in [3 p. 67-69], we deduce from 
the hypoelliptic estimates (43) and (44) that there exists some positive constants c 
and C such that 



(45) 
and 
(46) 



{z £ C : Re z > -1/2, Re z + 1 < c\z + l|^+r} n <j(q w (x, D x )) 



W^ixM-z)-^ <c\z + i\- 



2fc + l 



for all z £ C such that Re z > -1/2 and Re z + 1 < c\z + 1| 2fc o+ 1 . 

The property of exponential return to equilibrium may deduced from those resol- 
vent estimates. Indeed, Theorem |2.2| is a consequence of the abstract analysis led 
in [9] Theorem 6.1] or [15] (section 12, p. 754-756). Following more specifically the 
presentation given in [15] (section 12), we notice from (45) and (46) that under the 
assumptions of Theorem |2.2[ the quadratic operator q w (x,D x ) is a closed densely 
defined unbounded operator acting on L 2 (R n ) fulfilling the assumptions (12.1) and 
(12.2) in [H] (p.754) with 

1 
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Define for any t > 0, 



where 7 is a contour to the left of the spectrum that outside a compact set coincides 
with the curve 

Rez = hlmz\W«> +1 \ 
O 

with C > and 7 oriented in the direction of decreasing Im z. According to (45) and 
( 46 1 , this integral is convergent and defines a bounded operator depending smoothly 
on the parameter t > 0. It is then proved in [T5] that this integral operator is actually 
the semigroup e ~ tc T( x ^ D ^) generated by the accretive quadratic operator q w (x,D x ), 

2m .L, v 



Let < t < To with To the positive rate given by Theorem |2.1| Keeping on following 
the analysis led in |15] (p. 755) in the particular case when h = 1, we introduce two 
contours 7 and 7. Both contours are given by the curve 



Figure 2. 



Im z 




Rc z 



7ext\ 



Rez = I|i mz |i/(2feo+D ; 

in the region where Re z > b, with b — Re fXo +r. In the region, where Re z <b, the 
contour 7 is given by the equation Rc z = b; while 7 joins the two points 

b + iC 2k ° +1 b 2ko+1 and b - iC 2k ° +1 b 2ko+1 , 

further to the left so that 7 is entirely to the left of the spectrum of the quadratic 
operator q w (x,D x ), while 7 will only have the first eigenvalue in the bottom of the 
spectrum (1q given by Theorem 2.1 to its left, since < t < tq. Let 7j nt denote the 
vertical part of the contour 7 in the region where Re z = b and let 7 ex t denote the part 
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of 7 in the region where Re z > b. On the exterior piece 7 ox t, we have the resolvent 
estimate 

(47) \\(z-q w (x,D x ))- 1 \\ <0(\Imz\-W); 
whereas on the interior piece 7i n t, we have 

(48) \\(z-q w (x,D x )y 1 \\<0(l), 

since the compact contour 7; nt does not intersect the spectrum of the operator q w (x, D x ). 



According to Theorem 2.1 /io is an eigenvalue with algebraic multiplicity 1 for the 
operator q w (x,D x ). Arguing as in (12.13) in the article (15] . we may write that 



(49) 
(50) 



-^>^)= -L f e- tz (z-q w (x,D x ))- 1 dz 

= e-«*n Mo + _L f e- t *{z~q w (x 1 D x j)- 1 dz 1 



where stands for the rank-one spectral projection associated with the eigenvalue 
/io for the operator q w (x, D x ). Using the decomposition 7 = 7 int U 7 GXt together with 



(47) and (48), one can then estimate from the above the two parts of the integral 

1 



1 

2TTi 



as follows 



■■{z-q w {x,D x )) 1 dz = 



2m 



-tz 



Tint 



1 

2tt7 



(z-q w (x,D x )) X dz 

e- tz {z-q w {x,D x )Y 1 dz, 



1 

2tt7 



e~ tz (z - q w (x,D x )) L dz = L(L3) (e- bt ), 



7int 



2m 



e- tz {z-q w {x,D x )Y 1 d Z = O c{L 2 ) {l) / e -£» 2 *° +1 y~ W dy 
oxt Jc 2k o+ 1 b 2k o + 1 

/•+<» /I Is 

= O c(L2) (t- 2k °) e-*x 2k °~ 1 dx = O ciL > ) (l)(l + - + ... + -^)< 



-hi 



uniformly for t > 1. Since the accretive operator q w (x,D x ) defines a contraction 



semigroup, we deduce from (|49|), (j 50[) and the previous estimates that 
(51) 



W , + U C(L 2 )( 

uniformly for t > 0. Recalling that 6 = Re /io + t and considering now the operator 

Q = q w (x,D x ) - no, 



we directly deduce from (51) that there exists a positive constant C > such that 

Vt>0, || e -«2-n || £(L2) <Ce- rt , 

where Ilo is the rank-one spectral projection associated with the simple eigenvalue 
zero of the operator Q. 
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3.3. Proof of Theorem 2.3, Let q w (x,D x ) be a quadratic operator satisfying the 
assumptions of Theorem |2.1| Assume furthermore that this quadratic operator is 
real and satisfies ( 24 ) . As noticed in the discussion preceding the statement of Theo- 



rem [273} the lowest eigenvalue fiQ is necessarily real and the quadratic form a defining 
the ground state 



u (x) 



-a(x) 



G S(I 



q w (x,D x )u = HoU , 



is positive definite. It remains to check that the rank-one spectral projection Ho 
associated with the simple eigenvalue zero of the operator 

Q = q w (x,D x ) - /x , 



provided by Theorem 2.2 is orthogonal 

-2 



n M = ||wo|| L 2( K n)(u,wo)L2(M")«o, ueL (K n ), 



in order to deduce Theorem 2.3 from Theorem 2.2 This property is a direct conse- 
quence of the assumption 



(52) 



Ker(q w (x, D x ) - no) = Ker(q w (x, D x )* - Mo ) = Cu c 5(R n ). 



Indeed, notice from (6.52) in [STJ III. 6] that the adjoint of the spectral projection Ho 
associated with the eigenvalue zero for the operator Q is also a spectral projection for 
the operator Q* , 



(53) 



i 

2tTT 



1 

2Z7T 



{z-Q*)- 1 dz; 



with r being the mirror image of T with respect to the real axis and T a circular 
contour centered in the eigenvalue with a sufficiently small radius. Notice that 
in the integrals ( 53 ) both contours T and T are taken in the positive direction and 
that ITq is actually the spectral projection associated with the eigenvalue zero for the 
operator Q*. Since ±I and IIq are both projections, one may decompose the L 2 (K") 
space in two (possibly different) direct sums 

L 2 (M") = Ker(Ho) © Ran(II ) and L 2 (W 



Ker(n5) © Ran(n5) 



It follows from ( 52 ) that 



Ran(flo) = Ran(ITo) = Cuq. 



Indeed, Theorem 2.1 applies for both quadratic operators q w (x,D x ) and q w (x,D x )*\ 
and the eigenvalue zero has therefore algebraic multiplicity 1 for both operators Q 
and Q* . Furthermore, since 

Ker(IIo) = Ran(no)- 1 = (Ctto) 1 " and Ker(IIo) = Ran(n ) ± = (Cuo)^, 

we conclude that the spectral projection n is orthogonal 

n u = ||wo||7 L 2 2 (r „ ) (u,'«o)l2(r>»)Uo, u G L 2 (R n ). 

This ends the proof of Theorem |2.3| 



4. Some applications 
In this section, several applications of the singular space theory are discussed. 
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4.1. The Kramers- Fokker-Planck operator with quadratic potential. As pointed 
out in [16] , a noticeable example of quadratic operator with zero singular space is the 
Kramers- Fokker-Planck operator considered on the unweighted L 2 (M.^ v ) space 

,.2 



K 



-A, 



1 



+ v.d x - V x V(x).d v , (x,v) € 



with a quadratic potential 



V(x) 



-ax , a £ 



Indeed, this operator writes as 



with a symbol 



K = q w (x,v,D x ,D v )- 



q(x, v, V) = V 2 + \v 2 + i(v£ - axrj), 



which is a non-elliptic complex- valued quadratic form whose real part is non-negative. 
Simple algebraic computations show that the associated Hamilton map 

q(x,v,ti,rj) = a((x,v,£ > T]),F(x,v,t > T})), 

is 



F 



and that the singular space is equal to zero S = {0}. More precisely, the integer ko 
appearing in (22 1 is equal to 1: 

Ker(Re F) n Ker(Re F Im F)nl 4 = {0}. 

Now a computation shows that A £ C is an eigenvalue of 2F precisely when 



/ 













\ 




-\ai 








1 















\ai 




V 





1 

4 


-¥ 





/ 



A - 



A 



1 = 0, 



and we easily see that when a < 0, the eigenvalues A of 2F with Im A > are given 
by 



Ai 



When a > 0, we get the eigenvalues 



-i + i\J\ - 4a i + i\J\ - 4a 
^ ' A 2 = . 



i + Wl - 4a 
Ai — , A2 



i — iyl — 4a 



According to (18 1, the spectrum of the operator q w (x, v, D x , D v ) is given by 



•Ai 



Ai 

i 



A, 



kj £ N 



In particular, when a > 0, we observe that the lowest eigenvalue of the spectrum of 
the operator q w (x,v, D x , D v ) is 

_ 1 
Mo - 2' 

whereas when a < 0, this lowest eigenvalue is equal to 

VI -4a 
Mo = o • 
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The spectral gap tq > appearing in Theorem 2.1 is respectively equal to 



VI - 4a - 1 1 - -v/1 - 4a 1 
2 ' 2 ' 2' 

when 

a<0, 0<a<l/4, a > 1/4. 

Notice that the lowest eigenvalue is always real. This is consistent with the fact 
that the quadratic operator q w (x,v,D x ,D v ) is real. Further calculations allow to 
determine explicitly the ground state 

q w (x, V, D x , D v )u = (j, u . 

When a > 0, we obtain the usual Maxwellian 

u {x,v) = e -i M2 -T = e -i(£+v(*)) e S (M. 2 ), 

whereas when a < 0, the ground state is given by 



Notice further that the assumption 

Ker(q w (x, D x ) - Mo ) = Ker(<f (x, D x )* - Mo ) = Cii C S(M 2 ), 



holds true only if a > 0. Theorems |2.1| and |2.2| therefore apply whenever a / 0, 
whereas Theorem |2.3| only applies when a > 0. 



4.2. Chains of oscillators. This application comes from the series of works [2HH[3] 
and was studied in [T^] (Section 6). We shall begin by recalling the setting of the 
analysis |12j and explain how the results obtained there relate with the present work. 

This example is a model describing a chain of two oscillators coupled with two 
heat baths at each side. The particles are described by their respective position and 
velocity (xj,Uj) € M. 2d . For each oscillator j g {1,2}, the particles are submitted to 
an external force derived from a real- valued potential Vj(xj) and a coupling between 
the two oscillators derived from a real-valued potential V c (x2 — xi). We write V the 
full potential 

V(x) = V 1 (x 1 ) + V 2 {x 2 ) + V c (x 2 -x 1 ), x = ( Xl ,x 2 )eR 2d , 

y = (2/1,2/2) € R 2d the velocities and z = (zi,z 2 ) € R 2d the variables describing 
the state of the particles in each of the heat baths. In each bath, the particles are 
submitted to a coupling with the nearest oscillator, a friction force given by the friction 
coefficient 7 and a thermal diffusion at temperature Tj, j € {1,2}. We denote by Wj, 
j G {1,2}, two d-dimensional brownian motions with mean and variance 1; and 
write w — (wi,w 2 ). The system of equations describing this model reads as (see [6]), 

Vidt 

y 2 dt 

—d Xl V(x)dt + z\dt 
—d X2 V(x)dt + z 2 dt 
—jzidt + jxidt — ^/2jTidwi 
-jz 2 dt + jx 2 dt - v / 2jT 2 dw 2 . 



(54) 



dx\ 
dx 2 
dyi 
dy 2 
dz\ - 
dz 2 - 
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Following 12 and setting T\ — a\h/2, T 2 = a^h/2, the corresponding equation for 
the density of particles becomes 



(55) hd t f + l ai {-hd Zl )(hd Zl + — (zi - Xl ))f 
2 \ a.\ J 

+ ]:a 2 (-hd Z2 )(hd Z2 + - 
2 V a 2 



1 a 2 {-hd Z2 )(hd Z2 + ^-(z 2 - x 2 ))f + (yhd x - (V x V(x) - z)hd y )f 



Define 



and 



v 2 z 2 

<5>{x 1 y,z) = V{x)+ y - + - 



1 . 1 

M Q = —e 



with a > and C > 0. One may check that the function M ai with a = a\ — a 2 , is a 
Maxwellian of the process when the temperatures are the same. Considering now the 
general case when temperatures may be different, the function Ai a is used to define 
the weighted space L 2 (e~ ^ dxdydz) . In order to work in the flat L 2 space, we change 

1 /2 

the unknown by setting / = A4a u. The new equation for the unknown u reads as 

(56) hd t u+ ~ai(-hd Zt + -{z x -x^fhd^ + (— - -V«i - x 1 ))u+ 

2 \ a / V Qi a/ / 

^a^-hd^ + ^^-x^ (hd Z2 + ^ — ~) {z 2 -x 2 )^u+(yhd x -(V x V(x)-z)hdy^)u = 0. 

We consider the case when external potentials are quadratic. For simplicity, we may 
assume that /i = 1, 7 = 2, d = 1 and take 

(57) Vi(x 1 ) = ~axl, V 2 {x 2 ) = hjx 2 2 , V c (x 1 - x 2 ) = ^c{x 1 - x 2 ) 2 , 



with a, b, c € R. Equation ( |56| ) writes as 

dfti + D x )u - 2u = 0, X = (x, y, z) € K 6 , 

where q w (X, Dx) is the quadratic operator with symbol 

17 = aiC? + "zCl + Pi( z i ~ x if + fii{z 2 - x 2 ) 2 + i[25 1 C,i(zi - x\) + 28 2 C, 2 {z 2 - x 2 ) 
+ 2/1 £1 + 2/26 - m (( a + c ) x i ~ cx 2 - zx) - t]i ( - gei + (6 + c)x 2 - Z2)] , 

with 

a \ai a/ a \a 2 a/ a a 

the notations £, 77, £ standing respectively for the dual variables of x, y, z. The condi- 
tion 

1 

a > — max(ai, a 2 ), 

appearing in [12] (Section 6) exactly ensures that this quadratic symbol has a non- 
negative real part Re q > 0. Notice that only the case with identical temperatures 
a = oli = a 2 is discussed in 12J. More precisely, the authors mention that no 
Maxwellian is known in the case of different temperatures and that it prevents them 
from finding any supersymmetric structure. Here, we consider the general case with 
possibly different temperatures a% ^ a 2 for the quadratic potentials defined in (57 1 
and assume that 

a > — max(«i, a 2 ), a± > 0, a 2 > 0. 
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Direct algebraic computations using the definition of Hamilton maps and the explicit 
expressions 

Re q = aiCi + "2(2 + Pi( z i ~ ^l) 2 + £2(22 - x 2 ) 2 , 

Im q = 25iCi(^i - xx) + 26 2 ( 2 (z 2 - x 2 ) + + 2/26 

- 771 ((a + c)x\ - cx 2 — z\) — - cxi + (b + c)x 2 - z 2 ), 

show that (x, y, z, £, 77, £) = (Re F)(x, y, z, £, 77, () with 

it =0, x 2 = 0, yi =0, j/ 2 = 0,5i = aiCi, z 2 = CK2C2, 6 = - 

£2 = /?2(^2 - x 2 ), fji = 0, 772 = 0, Ci = ~Px{zi - xi)X 2 = -/3 2 (z 2 - x 2 ), 
and (x, y, z, £, 77, C) = (Im F)(a;, y, z, £, 77, C) with 



xi = -yi, x 2 = - 



1 V2, Vi = - \((a+c)xi- cx 2 - zi) , y 2 = -~(-cxi + (b+c)x 2 -z 2 ), 



z\ = &\{z\ - xi), z 2 = 5 2 {z 2 - x 2 ), |i = StCi + \i a + c ) 7 h - t,«72, 



I2 = S 2 ( 2 - t,c?7i + ^{b + c)rj 2 , 771 = -gfi, ?72 



1 



,£2, Ci 



-<5iCi - 2^1, 



C2 



It follows that 

Ker(Re F) n M 12 = {(x, y, z, £, ?),() 6R 12 :( = 0, ar = z}, 

Ker(Re F) n Ker(Re Flm F) n K 12 = {(x, y, z, £, 77, £) £ M 12 : y = rj = ( = 0, x = z} 
and 

Ker(Re F) n Ker(Re Flm F) n Ker(Re F(Im F) 2 ) n M 12 
= {y = £ = V — C — 0, x — z, (a + c — l)xi — cx 2 = 0, — cxi + (b + c — l)x 2 = 0}. 
When the condition 

(58) (a + c-l)(b + c-l)-c 2 ^0, 

holds, the singular space reduces to {0} after intersecting exactly ko + 1 kernels with 
here 0<fco = 2<ll. Notice that condition (58) corresponds exactly to the assump- 



tion V(x) — x 2 /2 is a Morse function required in [12] (Lemma 6.1) to ensure that 
the needed dynamical conditions hold. The other conditions d"Vj(x) = 0(1), when 
\a\ > 2, for j = 1,2, c; and |V x V(x) — x\ > 1/C when |x| > C, are also satisfied 



for q uadr atic pot entials fulfilling condition (58). When this condition holds, Theo- 
rems 2.1 and 2.2 applies. This shows in particular that the first eigenvalue]^] in the 
bottom of the spectrum /j,q has algebraic multiplicity one with an eigenspace 

(59) Kev(q w (X,D x ) - Mo ) = CM ai , aa , 
spanned by a ground state of exponential type 

(60) M ai , aa (x, y, z) = e-<*<»^ G 5(R 6 ), 



2. which may be computed explicitly by using the formula Jl9| 
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2:S 



where a is a complex-valued quadratic form on K 6 whose real part is positive definite. 
Notice that the operator 

1 , ,\ /„ ,2 1 



q w (X, D x )u = 2u + a x ( - d Zl + -( Zl - xx)) (d Zl + ( ) (zx - x 1 ))u+ 

\ a I \ ax a/ I 

• ( - d Z2 + \z 2 - x 2 )) (d Z2 + (A - (z 2 - x 2 ))u+ (yd x - (W x V(x) - z)d v ) >< 



a / \ ax a j 

(i 2 ( <>■_.. + -{Z2-X2)) (d Z2 + (—--){ 2 .>■■> 1 ! " - 1 !/<>, - 1 \ , > (•'•!- - i'/„ ) «. 

is real. The discussion preceding the statement of Theorem |2.3| then shows that the 
quadratic form a is positive definite. This proves the existence of a Maxwellian A4 ai .a 2 
in the general case when the temperatures may be different ol\ 7^ q 2 - It would be 
interesting to push further the computations in order to derive the exact expressions 
for the eigenvalue /iq, the spectral gap tq and the Maxwellian Ai ai ,a 2 - Another 



question of interest would be to investigate if Theorem 2.3 applies by checking the 



validity of condition (24). Those questions arc left for a future work. From now, we 



may only notice from ( 45 1 and ( 46 ) that the following resolvent estimate holds 
3c>0, {z G C : Re z > -1/2, Re z + 1 < c\z + n a(q w {X, D x )) = 0, 

30 0, \\(q w (X,D x )-z)- 1 \\ <C|z + l|-*, 

for all z G C such that Re z > —1/2 and Re z + 1 < c\z + since the integer fco is 
here equal to 2. 

4.3. The generalized Langevin equation. The spectrum of the Fokker-Planck 
operator associated to a linear stochastic differential equation in M d has been calcu- 
lated in [23], see also [23]. In this section we consider a particular system of linear 
stochastic differential equations that is obtained as a finite-dimensional Markovian 
approximation of the non-Markovian generalized Langevin equation (GLE) in R d , 

(61) x = -V x V(x)- f 7(i- s)x(s)ds + F(t), 



Jo 

where V(x) is a smooth confining potential and F(t) a mean zero stationary Gaussian 
process with autocorrelation function 7(f), in accordance to the fluctuation-dissipation 
theorem 

(F(t)®F(s)) = p- 1 j(t-s)L 
Here j3 > stands for the inverse temperature and I for the identity matrix. The GLE 
equation together with the fluctuation-dissipation theorem appear in various settings 
such as surface diffusion [T] and polymer dynamics |32j . This equation also serves 
as one of the standard models of non-equilibrium statistical mechanics describing the 
dynamics of a 'small' Hamiltonian system (the distinguished particle) coupled to one 
or more heat baths which are modelled as linear wave equations with initial conditions 
distributed according to appropriate Gibbs measures [59] . For further details, see [2"5] 
and the references therein. 

For simplicity, we will consider the case d = 1. As in |25j . we study the case when 
the memory kernel is a sum of exponentials 

m 

7 m (t) = E A i e " Qi|t| ) 



with aj > 0, Xj ^ 0, j — l,...,m. Equation (611 is then equivalent to the dynamics 
of the following system 

dx = ydt, 

(62) { dy = -V x V{x)dt + Y™ =x \ j z j dt, 

dzj — — (Xjy + o-jZj) dt + yj2ctj(3~ 1 dWj, j = 1, m, 
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with (x,y,z) £ E m+2 , z — (zi,...,z m ) 6 M m and {Wj)j=i... m being independent 
standard Brownian motions. 

We consider the specific case when V is a non-degenerate quadratic potential 

V(x) = 



with u =/= 0. The density of the invariant measure associated to the dynamics (62) is 
then given by 

with Zp being the normalization constant and | • | the Euclidean norm on M m . This 
process is Markovian whose generator 

m m m m 

(63) C = yd x - V'(x)d y + (E ^3) 9 j/ - 2/ E *jd Zj - E a i z i d H + E a o^ ld % > 
is an hypoelliptic operator acting on L 2 = L 2 (pdxdydz) , with L 2 -adjoint 

m m m m 

£* = -ya, + r'< - ( £ A^) d v + y E W„ + E (*i0 + E : '' 2 . 

3=1 1 = 1 3=1 3=1 

acting on By using the transformation 

H = p- 1/2 r( P 1/2 -), 

in order to work with an operator acting on the unweighted L 2 (1R™+ 2 ) space 

m mm ^ 1 m 

7* = - y a, +^xd y -(E A i^)a M +yE A A - E "i (- ^"X- + rl) + 2 E ^ 

3=1 3=1 3=1 3=1 

and setting 

(64) H = _h + -E«3, 

3=1 

we notice that 

m 

H = Q w (x,y,z,D xl D y ,D z ) = yd x - u 2 xd y + fE A 3%)^ 

3=1 

m 771 

3=1 3=1 

is a quadratic operator, whose symbol is given by 



L 2 

A' 



q{x,y,z,^,r],C) = i(y£ - u 2 xri) + i^X^Zji] - yQ) +^2a 3 \f3 1 C 2 +P 



z 2 



This quadratic symbol has a non-negative real part Re q > and a direct computation 
shows that its Hamilton map 

q(x,y,z,C,v,0 = (^({x,y,z,^,r],Q,F(x,y,z,^,r],C)), 

is given by (x,y,z, £,f}X) = F(x,y,z,£,ri,(), with 

1 11 »"\ 1 1 

x = ^iy, 5 = ~-^iw 2 x + ^i^XjZj, Zj = --iXjy + ajp~ x Q, £ = -iw 2 ij, 

3=1 
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2t, 



11™ 1 1 

3=1 



(XjZj 



It directly follows that 

Ker(Re F) n R 2 ( m+2 > = {(x, y, z, £, ry, C) € M 2 (" l + 2 ) : z = ( = 0}, 

Ker(Re F) n Ker(Re Flm F) n ]R 2 ( m+2 ) 

= {(*.!/, *>&»7,0 ^ 2K2) :y = z = V = ( = 0} 7 

Ker(Re F) n Ker(Re Flm F) n Ker(Re F(Im F) 2 ) n M 2 (" l + 2 ) 

= {(a:, y, z, v , Q & K 2(m+2) : x = y = z = £ = V = C = 0}. 

This proves that the singular space associated to the symbol q reduces to {0} after 
intersecting exactly ko + 1 kernels with here < fco = 2 < 2(m + 2) — 1. Furthermore, 
notice that the quadratic operator % — q w (x, y, z, D Xl D y , D z ) is real. Setting 

u (x,y,z) = Z X J 2 p{x,y,z) x l 2 = e -f(^) +3 # + ^) = e"! (™ 2 * 2 +^+M 2 ) e S (R™+ 2 ) 
and 



Mo 



^ 771 

3=1 



an > 0, 



a direct computation shows 



Hu = H*u = HqUq. 



On the other hand, notice that 

Re(Hu,u) L 2 =Re(H*u,u) L 2 = ^ a, {\\p~ 1/2 D Zj 

m 

>2^« j ||/3- 1/2 ^ 3 «I|l 



3=1 



3=1 



u\\ 2 T , + 



L 2 



> ^2aj\\D Zj u\\L» \\zju\\ L 2 > no\\u\\ 2 L2 , 



3=1 



since 



m m ^ m 



3=1 



3=1 



3=1 



and 



3=1 3=1 

Here the notation [D z .,iZj] stands for the commutator of the two operators D z . = 
i~ 1 d Zj and izj. This implies that /jq is necessarily the lowest eigenvalue in the spec- 
trum of the quadratic operator % and that it is the associated ground state appearing 
in Theorem |2.1| We notice further that all the assumptions of Theorems |2.1| |2.2| and 



2.3 hold true. Thus, these three theorems apply. This proves the property of exponen- 



tial return to equilibrium for this finite-dimensional Markovian approximation of the 
non-Markovian generalized Langevin equation with quadratic potential. Here, Theo- 
rem [2T] provides an exact formula for the rate of convergence that may be computed 
explicitly 

To = 2 min Im A. 

Xecr(F) 
Im A>0 
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For simplicity, consider the specific case when m — 1 and denote A and a for the 
parameters Ai and ct\. An explicit computation of the characteristic polynomial of F 
shows that 

P(X) = det(F-XI) = X 6 + ha 2 -2X 2 -2w 2 )X 4 + ^((X 2 +w 2 ) 2 -2w 2 a 2 )X 2 + ^-. 

4 16 64 

This is a polynomial equation of the third degree in the variable Y = X 2 which may 

be solved explicitly. Refraining from using Cardan formulas for giving an explicit, 

but rather complicated, formula for the rate of convergence r , we may instead try to 

compute numerically this rate in order to study its dependence on the parameters a 

and A. A limited attempt in this direction is provided by the following calculations 

where the rate tq is numerically computed for different values of the parameters A 

and a in the case when w = 1. Notice in particular that the spectral gap seems to 

be maximized along a curve of the type a = 7A 2 , with 7 > a positive constant. 

On the second figure, the spectral gap is computed as a function of the parameter 

A while keeping the friction coefficient 7 = ^- fixed. We observe that the spectral 

gap becomes almost constant for a large range of value s of A. This is in accordance 

with the numerical observations reported in [2]. It would be interesting to maximize 

the spectral gap, i.e. to optimize the rate of convergence to equilibrium, by choosing 

appropriately the parameters ctj and Xj in ( |62[ ) . This question, together with possible 

applications to Markov Chain Monte Carlo (MCMC) techniques will be studied in 

future work. 



Figure 3. Spectral gap tq as a function of a and A. 
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i 2 

Figure 4. Spectral gap t as a function of A with 7 = — fixed. 
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